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ABSTRACT 



This thesis presents a method for estimating the distribution of 
stress components around the outer circumference of a cross section of 
a pipe, from the readings of strain gage elements arbitrarily positioned 
and oriented around this circumference. A least-squares procedure is 
used to obtain best estimates of the coefficients of Fourier expansions 
describing such distributions. A digital computer program was developed 
for applying the method. Data for testing the method and program were 
generated by a computer program ufcing the best available theory of stress 
analysis in pipes. Methods for adding random errors to the data were 
adapted and used for closer simulation of actual situations. 

A second problem treated in this thesis is that of inferring the 
loading acting through the cross section of a straight pipe of concen- 
tric bore from known stresses at points on the external surface of the 
cross section which is presumed to be distant from stress concentra- 
tions. It is shown that this inference can be made from stresses at 
only three points of the cross section. A digital computer program was 
developed to do this. 
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1 . 



Introduction. 



In a thesis presented to the Naval Postgraduate School in May, 1966, 
LT Joseph W. Koch, Jr., USN, did pioneer work on the problem of obtain- 
ing optimal inferences of stresses in submarine main sea water piping, 
from strain gage readings taken during dockside hydrostatic tests and 
during submerged operations. This work provided a sequential procedure 
for determining such stresses, but it also indicated the need for many 
related investigations . 

Two of these are treated in the present thesis. The first is that 
of inferring, in an optimal way, the distributions of stresses around 
the outer circumference of a given cross section of pipe from observa- 
tions of strain obtained from an arbitrary distribution of strain gage 
elements around that circumference. There are several reasons for being 
interested in this problem. One reason is that even though the gage 
elements may be arranged in special patterns intended to simplify deter- 
mining stresses (for example three-gage rosette arrangements), failure 
of any of the elements might make it difficult to make effective use of 
valid information obtained from closely related gages. Another reason 
is the ability to infer stresses from arbitrarily arranged gages, which 
is a logical step in investigating the most effective and most economi- 
cal arrangements of gage elements. 

Briefly, this problem is solved herein by constructing Fourier 
series representing the desired stress components. Theoretical consider- 
ations are used to evaluate the coefficients of the series in an optimal 
way, minimizing the sum of the squares of the residuals. A digital com- 
puter program was developed to support the theory. 
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The second problem treated here is that of inferring the force and 
moment components applied through a cross section of pipe and the inter- 
nal pressure acting inside the pipe. It is shown that these values may 
be inferred from a knowledge of the state of stress at only three abritrar- 
ily chosen points at the external circumference of the cross section, and 
a digital computer program for making the corresponding evaluations was 
developed. 

The solution of this problem assumes that the stress distribution is 
that given by what may be regarded as the best available current theory 
for the state of stress in a section of uniform pipe having concentric 
bore and remote from stress concentrations and other disturbances in 
the stress pattern. 

A digital computer program to generate numerical data from the best 
available theory was developed. This data was used to test and verify 
the theory and the computer programs. Suggestions of how to extend the 
present work are given in Section 4. 

The author would like to express his sincere appreciation to Dr. 

John E. Brock, Professor of Mechanical Engineering, Naval Postgraduate 
School, for his continued patience, efforts and most capable guidance 
while acting as faculty advisor in this work. 



6 



2. Least-Squares Estimation of Stress Distributions at the external 
Circumference of a Pipe by Means of Fourier Expansions, Using 
Strain Gage Data. 

The usual method of estimating the stress distribution at 
the surface of a body is by taking readings of strain gages grouped in 
rosettes, to determine local principal stresses. Unless the strain 
rosettes have a redundance of gage elements, when one or more of them 
fail, it is difficult to make effective use of the information provided 
by the elements which do not fail. 

The work that follows presents a method to estimate the stress 
distribution around an external circumference of a cross section of pipe, 
along a sequence of strain gage elements, which may or may not be group- 
ed in rosettes. If one or more gage elements fail, the information 
from the others is still useful. 

It is assumed that strain gage elements are arbitrarily placed 




AA view 




(a) (b) 

Fig. 2.1. 

and arbitrarly oriented around the external circumference of a cross 



section of pipe so as to measure the strain produced by application of 



loading. Assume that there are N such elements as indicated in Fig. 

2 . 1 . 

Define two orthogonal directions at the cross section, one normal 
to it, called longitudinal and represented by the letter L, and the 
other tangential to the external circumference, called tangential and 
represented by T; Fig. 2.1.b. Also define the angle 9 as the angle be- 
tween the vertical line that passes downward through the center of the 
cross section, and the radial line that passes through any specified 
point of the external circumference, measured as indicated in Fig. 2.1.b 

Measure the angular orientation of each gage element in the counter 
clockwise direction from the tangential axis T, and call this angle 0. 

A set of parameter-pairs can then be formed 



(0. , 0.) i “ 1, 2, ... , N 

which describe the position and orientation of each element, Fig. 2.2. 




Fig. 2.2. 

At each one of the gage locations, there are tangential, longi- 
tudinal, and shear strains, (€•).,( e ) 

, T i L l ~ 

and ( €, ) . respectively. The strains 

1 L i 

€ and € along axis y and z, rotated 
y z 



the angle 0 from the tangential direction 
Fig. 2.3., are related to the local 




Fig. 2.3. 
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strains € , € and € by the strain tensor in the following way 

i Li 1 Li 



€ %« 1 

y yz 


i 


cos 0 sin 0 




** 

e \ e 

T ^ TL 




cos 0 -sin 0 


h 6 « 

yz zj 




-sin 0 cos 0 

, i 




^ TL L 

V 1 




sin 0 cos 0 

k i 



The corresponding relations for the strains ( € and ( * z )^ 
along and transverse to the ith strain gage axis, are then 

( % } i “ (« T ) i sin 2 0 i + (e^^os 2 ^. + (« TL ) iS in 0 iC os ^ (2.1.1) 

( * 2 > i = ( e T ) i cos‘ i 0 i + ( € L ) iS in 2 0 - ( « TL ) i sin0 i cos (2.1.2) 

At the external surface, strains are related to the corresponding 
stresses as follows 



* T ■ < S T - V S L> 


/E 


(2.2.1) 


* L ‘ (S L - V S T > 


/ E 


(2.2.2) 


6 T L " S TL ' G 




(2.2.3) 



where s is stress 

v is Poisson's ratio 
E is Young's modulus of elasticity 
G is the shear modulus of elasticity 
Substituting these expressions for the strains in the relations 
(2.1) for the strains along and transverse to the ith gage element and 
simplifying 

(€ y ) i = (s T ).[l - (v + l)sin 2 0 i ]/E + (s L ) i [l - (v + l)cos 2 0 l ]/E 



+(s TL ).sin 2 0./(2G) 



(2.3.1) 
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(2.3.2) 



(€ z > i = (s T )[l - (v + l)cos 2 0 i ]/E + (s L ) 1 [l - (v + l)sin 2 0 i ]/E 

- (s TL ).sin2 0./(2G) 

A strain gage element is sensitive not only to the strain along 
its own axis but also, to a small extent, to the strain at right angles 
to this axis. The transverse sensitivity is usually hard to obtain. It 
is rarely available from the manufacturer and may have to be inferred from 
an experiment. Nevertheless, it may be of importance and it is included 
in the theoretical development that follows. 

For a particular gage element properly mounted on a given definite 
material, the relative change in electrical resistance, which is sensed 
by the strain indicator (or recorder), is given by a formula of the sort 

Ar = F € + F€ (2.4) 

— y y 2 * 

where R is the electrical resistance of the element, € and £ are 

y z 

the strains along and perpendicular to its axis, and F^ and F z are co- 
efficients which depend on the characteristics of the gage and the pro- 
perties of the material upon which it is mounted. 

The ratio 

k = F / F (2.5) 

z y 

is called the transverse sensitivity (factor) of the gage, and the 
quantity 

F' = (1 - k 2 ) F 

y 

is called the gage factor (which should be set on the sensor instrument). 
F' is the ratio of unit change in electric resistance to the quantity e, 
(indicated strain) sensed by the instrumentation. 
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Thus 



) / F' 



- ( F y e y + F z € z )/[(l - k 2 )F y ] 

= ( € + k€ z )/(l - k 2 ) (2.7) 

This is what the instrumentation should indicate if there were no 
errors. Clearly for a gage with no transverse sensitivity, k = 0, and 

e = € . 

y 

After substituting the expressions (2.3) for ( € y )^ and ( €. z )^ in 
equation (2.7) and simplifying, the theoretical indicated strain e^ 
sensed along the ith strain gage element is given by 

®i = ( S T^i^ X ~ kv) - (1 - k)(l + v)sin 2 0 i ]/[E(l - k 2 )] + 

+ (s L ).[(l - kv) - (1 - k)(l + v)cos 2 0 i ]/[E(l - k 2 )2 + 

+ ( S TL ) t < 1 - k)sin20 i /[2G(l - k 2 )] (2.8) 



Let the coefficients of the stresses in the above expression be 
represented by 

A(0 i ) - A t - t(l - kv) - (1 - k)(l -fv)sin 2 0.]/[E(l - k 2 )] (2.9.1) 
B(0 £ ) = B. = [(1 - kv) - (1 - k)(l + v)cos 2 0 i ]/tE(l - k 2 )] 

( 2 . 9 . 2 ) 

C(0 t ) - C ± = (1 - k)sin20 i /[2G(l - k 2 ) ] (2.9.3) 

Then equation (2.8) can be put in the form 



e i A i (8 T ) i + B i (9 L^i + C i (i Wi 



( 2 . 10 ) 



Any one of the stresses (s_,) . , (s T ), and (s„ T ) has continuous and 

lX LX XL 

finite values around the external circumference of the pipe cross 
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section, and can be represented as in Fig. 2.4. 

The functional relation between stress and the angular position 
0 can be expanded in an infinite Fourier series^ - 



oO 



(s T )i “ 


'HiO 


(8 TA>n COS 


n© i + 


DU 

Si 

>1-0 


^‘Wn 


sin n0^ 


(2.11.1) 


< a l?i = 


oo 

>i:o 




n0 £ + 


Cm* 

o 


^LB^n 


sin n0. 

l 


(2.11.2) 


* S TL^i = 


oo 

s: 

hs Q 


( W« CM 


n0^ + 


H- o 


(S TLB^ 


sin n0. 
n l 


(2.11.3) 


where 
^ S TA^n ’ 


^ S TB^n ’ 


‘Vn • 


(S LB^n 


, (s 


SHA^n 


’ and < s SHB>n are 



as yet unknown numerical coefficients. 

Sufficiently accurate values can be obtained if the expansion is 
carried to a finite number of terms, say P terms, in which case the 
index n would reach a maximum value of P/2. Call this value M. Then 

‘Vi =■ ■§; ‘Vn™ 1 " e i + fr. <s TB ) „ sln n6 i <2 - 12 - 1) 



The conventional expansion of a Fourier series is a constant term 
followed by a summation from 1 to oo . The same is obtained by eliminat- 
ing the constant term and making summations from 0 to oo . This latter 
form will be used in this work to facilitate mathematical manipulations 
and later programming for a digital computer. 
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( 2 . 12 . 2 ) 



(s L>i 



M H 

(s TA ) cos n9. + «2L* (s to ) sin 9. 

Hso LA n i ' LB n i 

M /¥ 



<%>! » <Wn C0S n9 i + f? c ^TUSV 1 " n9 l < 2 ' 12 - 3 > 

Note that this expansion contains 6(M + 1) Fourier coefficients. 

If these expressions for the stresses are substituted into equation 
(2.10) for the strain e^, theoretically sensed by the ith strain gage 
element. 



e i *** 


S 

A . [ ^ 

i >?zo 


< s TA>n COS 


n9. 

l 


+ s 




+ V 


M 

2-1 • « • 


+ 


A/ 

z: . 



. . . ] +■ 



(2.13) 



H 



M 



+ C [ JEi ... + ZL (s ) sin n9 ] 
i * * o n i 



Call e* the strain actually measured by the ith strain gage ele- 
ment, and E. the difference between that measured strain and e., the 
l l 

strain that theoretically it should have measured, so that this error 
is expressed as 



E = e. - e* 
ill 



(2.14) 



the sum of 


the 


squares 


of the errors. 


that 






* 








R = 

A/ 


is# 


(E f ) 2 




= 


xc # 


< e i - 


e*) 2 




li 

R* 


(e.) 2 


- 2 


N 

e.e* + 


st 


<zt 


i 




/ i i 


VI 1 



; e *) (2.i5) 



To minimize R, partial derivatives respect to each of the 6(M + 1) 
Fourier coefficients are taken and set equal to zero: 



^ R ' 0 
^ D k 

where D fc = < s TA^k 1 (s TB>k • <S LA>k 1 <*LB>k 1 <S TLA>k • < s TUS>k in tur " 
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for k = 0, 1, 2, ...,M 



Taking derivatives of R from expression (2.15) 



R 



0 = 2 



C> 



Simplifying this 



*1 e 
<9i i 'b 



6 i 



V 










V 


S3 









6 i 



- 2 S.* e* 

i 



2 e i 

e ,. 



2_!j 
■J D. 



Taking derivatives of e^ from expression (2.13) 



* e . 



i ( s TA ) k 



= A^cos kQ^ for k = 0, 1, 2 M 



> e i 

* < S TB>k 



A sin k9. 
i i 



^ 6 i 
^ ^LA^k 

3 6 i 

^ S LB^k 



B.cos k9. 

l l 



B, sin k9. 
i i 



* 6 i 

^ ^ S TT.A^ 



TLA k 
i 



^ e 



C^cos k9^^ 



^ ^®TLB^k 



C.sin k9. 
l i 



(2.16) 

(2.17) 

(2.18.1) 

(2.18.2) 

(2.18.3) 

(2.18.4) 

(2.18.5) 

(2.18.6) 



Note that there are 6(M + 1) derivatives of this form. Substitut- 
ing in equation (2.17) the expression for e^ from equation (2.13) 



N r r 2 j —i 

1 L ^tiA 8 * 11 nG i + •♦•J + ••• + c i £ ••• 

^ S TLB^n Sin n9 i Ij ^D k ^ & t $ D 



... + 



i 

k 

(2.19) 
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which is the same as 



Af *Y 

\ fe A i 000 * 



*•50 




-J + 



Af rs* 

L Hfi***** 




e* 2i^ 

4 ^ Or 



( 2 . 20 ) 



for k = 0, 1, 2, ... , M 

When the index n in the Fourier expansions is equal to zero, the 

terms containing sin(n9 ) are zero. This makes the coefficients (s—,).., 

i ru u 

(s__) n and (s„, T _) n also zero, reducing the Fourier coefficients to (6M + 
LB U 1 Ld U ^ 

9 e . 

3) and the same number of derivatives of the form ■ ■ . There will 

k 



be (6M + 3) equations in 'each of which appear all the Fourier coeffici- 
ents. This makes a set of (6 m + 3) equations with the same number of 
unknowns, (the Fourier coefficients), which can be solved. 



Once the Fourier coefficients are solved for, they can be substitut- 
ed in expressions (2.7) for the stresses. The stress distributions s^,, 
s T and s are defined then as functions of the angle 0 over the com- 
plete external circumference of the cross section. 

A digital computer program was developed to apply the theory of the 
preceding text. It was tested using a mathematical model of a canti- 
levered pipe. The program is called TERESITA and is presented in Ap- 
pendix B of this thesis. 

It was assumed that strain gages having a zero transverse sensiti- 
vity factor were placed at the external circumference of cross section 
No. 1, see Fig. A. 1. in Appendix A. Four arrangements of gage elements 
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and several numbers of Fourier coefficients were used. 



The first arrangement consisted of twelve strain gage rosettes 
equally spaced 30° from each other around the cross section. Each rosette 
has three gage elements oriented at 0? 45°and 90°from the positive tangen 
tial axis. Runs were made with 11, 9, 7, 5 and 3 Fourier coefficients 
for each stress component. 

The second arrangement consisted of twelve rosettes spaced as above 
but with elements oriented at 0°, 60°, and 120°. Runs were made for the 
same number of Fourier coefficients. 

The third arrangement consisted of eight strain rosettes equally 
spaced 45° from each other. The rosettes had four gage elements orient- 
ed at 0°, 45°, 90° and 135° from the positive tangential axis. Runs were 
made for the same number of Fourier coefficients. 

As input data for testing the program, stress distributions pro- 
duced by the loading were evaluated at the external circumference of 
cross section No. 1. From the stresses, the theoretical strain that the 
gage element would sense, was computed. The stress distributions re- 
covered by the program agreed in all cases with the input information. 
Results using eleven Fourier coefficients (for each stress component s^, 
s T and s ) are almost identical with those using only three coefficients 

Additional runs were made using the second arrangement of gages, 
but eliminating arbitrarily the readings of some of them. This is what 
would have to be done in an actual case if some gage elements were known 
to have failed. Results are satisfactorily close to the preceding ones. 
Other runs were made using a strain gage transverse sensitivity factor 
of 1%. The program recovered likewise the input information. 

The fourth gage arrangement had one gage element every 10° around 
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the circumference with different orientation for each gage. Runs were 
made using the same number of coefficients. 

All the runs described above were repeated introducing in the gage 
element readings small random errors so to simulate actual situations. 
Calculated loadings differed from the input loadings to a degree which 
depended in each case on the amount of error introduced. 

The program and some sample numerical results of the runs made are 
presented in Appendix B of this thesis. 

A more detailed discussion of the effects of changing the number 
of Fourier coefficients and gage element arrangement is presented in 
Appendix E. 
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3. Inference of the Loading Acting at a Cross Section of a Pipe of 
Concentric Bore from Stresses at the External Circumference of 
that Cross Section. 

If stresses are known at a sufficient number of points around 
the external circumference of the cross section of a pipe of concentric 
bore at a distance from stress risers and concentrated loads, the load- 
ing acting at the cross section can be determined from theory. 

It is intended to analyse the effects that such loads produce 
at a general point P under the conditions specified above. From this 
analysis, general conclusions will be reached which lead to a mathe- 
matical procedure to determine the loads. 

Consider a cross section of pipe as specified above, Fig. 3.1. 
Choose reference axes x, y and z with the center point 0 as origin. Let 




Fig. 3.1. 

I, J and Tc be unit vectors along the positive x, y and z axis so to form 

a right handed triad. Let e be a unit vector along the outwards direc- 

tion of the radius that passes through the point P, and e^ another unit 

vector so that e , e and k form a right handed triad. The angle 0 is 
R T 
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measured from i counterclockwise to e , Consider the differential area 

R 

surrounding the point P. Assume that tensile stresses s and s (along 

L Ju 

the directions of e and 1c respectively) and shear stress s , are pro- 

1 TJj 

duced by internal pressure p and by a general force F and a general mom- 
ent M acting on the cross section at the center point 0. 

Let the components of F and M in the e , e and k triad be (F , F , 

Hi r r 

F ) and (1L, M ) respectively, and in the i, j, k triad, (F , F , F ) 
z K l z x y z 

and (M , M , M ). 
x y z 

From strength of materials considerations reviewed in reference 5, 
the stresses produced by p, F and M at the external point P are: 



o 2 , , .2 2v 

s T = 2pa / ( b - a) 



--Hjb/I+F^/A 



s tl - M z b / (2 I) + D F T 



(3.1) 

(3.2) 

(3.3) 



where a = internal radius of the cross section 

b = external radius of the cross section 
A = cross sectional area 

I = moment of inertia of the cross sectional 
area about a diameter 

D - [2 n a 2 + (2n - 2)b 2 ] / [4I( 1 + v)] 

v = Poisson's ratio 

n » v + 3/2 

Making use of the two dimensional coordinate transformation tensor 



“r 









cos 0 sin 9 

-sin 9 cos 9 



M 



M 



M can be expressed as 
T 



M = - M sin 9 + M cos 9 
T x y 



(3.4) 
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Similarly, 



F m = -F sin 9 + F cos 9 (3.5) 

T x y 7 

The expression (3.2) for s can be transformed to 

L 

s. = (1/A)F + (b/l)sin 9 M - (b/I)cos 9 M (3.6) 

L z x y 

Multiplying through by (I/b) 

(I/b)s - (l/(Ab))F + (sin 9)M - (cos 9)M (3.7) 

L z x y 

Letting (I/b) = g, and (I/(Ab)) = h, both of which are cross 

sectional constants, substituting above and changing sides 

h F + sin 9 M - cos 9 M - g s_ (3.8) 

z x y l 

The preceding theory is applicable to a point on the exterior 
surface of a piece of straight, uniform, concentric pipe at a suf- 
ficiently great distance from its ends or other connections which might 
introduce local disturbances into the stress distribution. Presuming 
that these conditions are satisfied, we next turn to the matter of try- 
ing to infer the loading conditions at the cross section, i.e., F , F , 

x y 

F^, M^, M z and p, from a knowledge of the stresses at several points 

on the external surface of the cross section. 

For a general point P, s. and 9 are known and the only unknowns in 

Li 

equation (3.8) are M^ and M^. Taking three such points at the cross 
section, the three unknowns can be determined proceeding as follows. 

Assume that the angular positions for the points are 9^, 92 and 
9^ and the axial stresses are s^, S 2 and s^ respectively. A set of 
equations can then be formed 

h F + s9. M - c9, M - gs. (3.9.1) 

z lx 1 y e 1 

h F + s9„ M - c9„ M ■ gs 0 (3.9.2) 

z 2 x 2 y 2 
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h F + sO. M - c9_ M = gs. 
z 3 x 3 y ° 3 



(3.9.3) 



where the letters s and c with no subscripts are abreviations for sin 
and cos respectively. These abbreviations will be used throughout this 
section. 

The set of equations (3.9) is solved applying Cramer's rule 



Denominator = 



s9. 



-c9. 



Numerator of F_ = 



Numerator of M 



x 



l v 3 



Numerator of M 
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s9 2 


C0 2 


h 


s ®3 


c6 3 


L - 9 


2 > + S(0 2 


-v 


8S 1 


S6 1 


-c9 


gS 2 


s9 2 


- C9 2 


gs 3 


s9 3 


-c9 3 


^2 " 


e 3 ) + s 2 


s(0 3 - 
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gS l 


- C9 1 
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8 8 2 


' c9 2 
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gS 3 


' c9 3 


• s 2 ) 


+ c9 2^ S 1 


- S 3> 
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S9 1 


gS l 


h 


s©2 


gs 2 


h 


s9 3 


gS 3 


' S 2> 


+ s0 2 ( s i 


- » 3 ) 



3 v 1 2" 

(3.11) 



3 v 2 



(3.12) 



(3.13) 
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The values of and are found by dividing the respective 

numerator by the denominator. 

The remaining components of the vectors F and M appear in the 
expression (3.3) for s^j which by using equation (3.5) can be trans- 
formed to 

s mT = (b/2I)M - Dsin 9 F + Dcos 9 F (3.14) 

TL z x y ' 

Multiplying through by (1/D), letting (-1/D) = g*, (-b/2ID) = h*, 

both cross sectional constants, substituting in (3.14) and changing 

sides 

h* M + sin 9 F - cos 9 F a g*s_ T (3.15) 

z x y ° TL v ' 

This expression is of the same form as (3.8) which led to the 

evaluation of F z> and by using information from three points. In 

this case the unknowns are M , F and F which can be evaluated by 

z’ x y 

using the same three points and applying the same method. In this form 
the loads F and M are completely determined. 

Since s^, is produced by the internal pressure p alone, it can be 
solved for from equation (3.1) with the information from any one of 
the three points used to evaluate F and M 

p = s T (b 2 - a 2 )/(2a 2 ) (3.16) 

From the above development, it can be concluded that theoretically 
the loads F, M and p acting at the cross section can be determined from 
known stresses at three points of the external surface of the cross 
section. 

A digital computer program called TIBISAY was developed to apply 
the preceding treatment. The program was tested using a mathematical 
model of a cantilevered pipe as in Section 2 of this thesis. Stresses 
generated at cross section No. 1, at points 0°, 120° and 240° from the 
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vertical downwards direction were used. The recovered loads were almost 



identical with the input ones producing the stresses. Deviations are 
due to small round-off errors introduced by the computer. The program 
and results of the test are presented in Appendix C of this thesis. 
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4. 



Conclusions . 



The results obtained in the tests for applicability of the theories 
developed in this thesis indicate that the proposed methods can be used 
for analysing stresses in pipes. 

There are, however, many factors influencing the validity of the 
results. Some of them are related to the amount of error introduced 
when determining geometries and material properties. Others are due to 
proximity to concentrated loads and stress concentrations. Other fac- 
tors are directly related to the degree of redundance of strain gage 
data, the number of gages used and their arrangement around the cross 
section, and the number of Fourier coefficients to be used. Also, it may 
be that the least-squares method used does not produce the best results; 
perhaps, by minimizing some other function of error, better results could 
be obtained. 

First steps were taken to investigate further some of these factors. 
Appendix E of this thesis presents results of a survey made concerning 
the influence of the number of Fourier coefficients used in the accuracy 
of results. There are also preliminary results of the influence of the 
number of gage elements and their arrangement around the cross section. 
Due to time limitations, further investigations were not conducted. 

The digital computer programs included in this thesis can be used 
to investigate deeper into the field. Some refinements of the methods 
can be added. It is believed that the following points should not be 
neglected. : 

a. Other criteria for minimizing errors can be investigated. 

A discussion of such criteria is given in reference 7 of this thesis. 

It is likely that some of the methods described there may have applica- 
tion in the pipe stress problem considered herein. 
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b. One may investigate the question of the minimum number of 
strain gage elements necessary to obtain acceptable stress distributions, 
and their orientation around the cross section. Similarly, the optimum 
number of elements and preferred orientation could be investigated. 

c. While observing the results for a determined gage arrange- 
ment, it was noticed that good results were obtained using from 3 to 7 
Fourier coefficients, but more than that produced bad results; using 
again 5 coefficients and changing the gage arrangement, good results were 
obtained for some arrangements and bad results for others. Some of these 
results are presented in Appendix G of this thesis. Further study could 
be made of the optimal number of Fourier coefficients to be employed. 

d. A method for data rejection should be incorporated in the 
theory and computer programs. There are instances when it is obvious 
that a strain gage element has failed and the reading of that gage should 
be neglected. Usually bad gages are not easily identified and bad data 

is processed along with good data. There is the possibility of incorporat 
ing a data rejection procedure in the method developed in Section 2 of 
this thesis. One suggestion for such a procedure is given in the follow- 
ing. 

The estimated Fourier coefficients of the stress distributions are 
evaluated from a minimization of the sum of the squares of the errors, 
i.e., the differences between the gage element readings and what they 
should have read theoretically. We can assume that these errors are 
random having a statistical normal distribution with mean and standard 
deviation which could be computed. Establishing a level of significance 
for a deviation to be sufficiently apart from the mean, one can reject 
the data having greater deviations and recompute the stress distribu- 
tions with the remaining data. Similar screenings can be done with the 
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newly computed stress distributions until the remaining data have de- 
viations under the level of rejection. This theory could be readily 
used if the redundance of data were large. Usually this is not the case 
and one must avoid rejecting data which are not perfect but which are 
useful and necessary. If the data is known to have a normal curve with 
a very sharp peak, bad data is easily determined and low significance 
levels can be used. If the bell is fairly flat, the majority of the 
deviations have pronounced deviations from the mean and the level of 
rejection becomes difficult to establish. The writer is inclined to 
think that there are criteria that may satisfy all conditions. Once 
they are found, the theory could be applied to this problem and to 
similar ones. 
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APPENDIX A 



DATA GENERATION THEORY AND RELATED DIGITAL COMPUTER PROGRAMS 
The data used to test the applicability of the theories developed 
in this thesis was generated for a mathematical model of a cantilevered 
pipe of concentric bore subjected to internal pressure p and general 
force F and general moment M applied to the tip. Details of this pipe 
are shown in Fig. A.l. The following material properties, pipe dimensions, 
and loads were specified: 





— y 




Fig. A.l. 

E = 30,000 ksi 
v * 0.3 

3 

Specific weight = 0.283 lb/in 
External diameter = 12.0 in 
Internal diameter = 10.0 in 



F 


= F 


= F 


= 100.0 


lb 


X 


y 


Z 




M 


= M 


= M 


= 100.0 


lb -in 


X 


y 


z 





p = 100.0 p6i 

Subroutine CANTIL presented herein in this Appendix, was used 
to evaluate the tangential, longitudinal and shear stress distributions 
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at the external circumference of cross sections 1 and 2 of this pipe. 

The stress distributions were used in turn to supply the data required 
for the tests specified in other parts of this thesis. 

This subroutine is used by either of the (main) programs TERESITA 
or TIBISAY described in Appendices B and C. The arrangement of the card 
input is indicated in Fig. B.l. of this thesis. Description and specifi- 
cations for CANTIL follow this page. 
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DESCRIPTION OF SUBROUTINE CANTIL 



Subroutine CANTIL evaluates the state of stress at any point of 
a cantilevered straight pipe of concentric bore, with internal pressure 
and loaded at the tip by general forces and general moments. The program 
is capable of analysing from 1 to 360 points in each of 1 to 5 cross 
sections of the pipe. The points have to be at the same radial distance 
from the center of the pipe, and separated from one another by, angular 
intervals which are submultiples of 360°, from 1° to 360°. 

The data deck, see Fig. B.l, includes the following cards: 



Card No. 


1. Format 3F20.8 


a. 


External diameter of the pipe, inches. 


b. 


Internal diameter of the pipe, inches. 


c. 


Pipe length, inches. 


Card No. 


2. Format 3E20.8, 2F10.8 


a. 


Young's modulus of elasticity, psi. 


b. 


Shear modulus of elasticity, psi. This can be 
omitted in which case it will be computed by the 
program from E/2( 1 + v) . 


c. 


Coefficient of thermal expansion. This can be omitted. 
It is not actually used in any present application. 


d. 


Poisson's ratio. 


e. 


Specific weight of the material, pounds per cubic 
inch. 


Card No. 


3. Format 215, 5F10.5 


a. 


Number of points to be analysed in each cross section. 


b. 


Number of cross sections to be analysed in the pipe. 


c . 


The distances from the tip of the pipe to the cross 
sections to be analysed, inches. 
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Card No, 4. Format F80.2 

a. Radial distance of the points to be analysed, inches. 
NOTE: for applications to other parts of this thesis, 
evaluations are made only at the outside sur- 
face, but the program is capable of making 
evaluations at interior points. 

Card No. 5. Format F80.2 

a. Internal pressure, psi. 

Card No. 6. Format 3F20.8 

a. Components of the force applied at the tip of the 
pipe, vertical, sidewards and outwards from the 
tip, so as to form a right handed set of orthogonal 
axes, lbs. 

Card No. 7. Format 3F20.8 

a. Components of the moment applied at the tip of the pipe, 
in the same arrangement as for the forces, lb-inches. 

Card No. 8. Format F80.2 

a. This is a sentinel card with one number on it. A 

negative number indicates that the job is finished; 
a zero or any positive number indicates that the 
program will be rerun for another set of radial 
distance and loads, in which case, cards similar to 
No. 4 through 8 must follow this one. 

OUTPUTS (all arguments are outputs) 

RO. External radius of the pipe. 

RI. Internal radius of the pipe. 

XE. Young’s modulus of elasticity. 
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G. 


Shear modulus of elasticity. 


XNU. 


Poisson’s ratio. 


SPECW. 


Specific weight of the material. 


D. 


A constant evaluated when analysing shear stresses pro- 



duced by shearing forces. 



XLl. 


A one dimensional array containing the distances from 



the tip to each cross section considered. 



NPOINT. 


Number of points analysed in each cross section. 


AREA. 


Cross sectional area of the pipe. 


ALL1. 


A three dimensional array intended to store the stresses 



at any of the points, in case the points are at the 





external surface of the pipe; otherwise, it has no mean- 




ing. 



Other subroutines which are needed and listings for which are given in 
this Appendix are: 



TORS: 


used by CANTIL to compute shear stresses produced by 



torsional moments. 



BEND: 


used by CANTIL to compute bending stresses. 


AXIAL: 


used by CANTIL to compute stresses due to axial forces. 


SHEAR: 


used by CANTIL to compute shear stresses produced by 



shear forces. 



CROSS : 


computes the cross product of two vectors. 


DOT: 


computes the dot product of two vectors. 


TIMES: 


computes the product of a scalar and a vector. 


PLUS: 


computes the sum of two vectors. 


MINUS: 


computes the difference between two vectors. 


TRANSF: 


computes the components of a vector in another angularly 




displaced set of axes. 
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METHOD. 



The vectorial approach discussed in reference 5 of this thesis is 
used to compute stresses produced at points of the cross section by 
internal pressure and forces and moments acting at the centroid of that 
cross section. 




Fig. A-2. 

Those forces and moments, in turn, are computed by the classical vectori- 
al analysis of the statics of a beam. If the force F and the moment M 
are applied at one end of a beam, Fig. A. 2. a., the force F^ and moment 

M acting at the centroid of the cross section nn, Fig. A.2.C., are 
n 

found by considering the equilibrium of the segment of beam from the 
tip to the cross section. Summing forces and taking moments about the 
centroid of the cross section 
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F + W + F n = 0 (A. 1) 

M n + M + LxF + %LxW = 0 (A. 2) 

where L is the vector distance from the centroid of the cross section 

to the point of application of the loads at the tip, and W is the force 

due to the weight of the segment of beam. 

Solving for F and M 
n n 

F = - (F + W) (A. 3) 

M n = - (M + LxF + %LxW) (A. 4) 

The forces and moments acting at the cross section under considera- 
tion on the remaining portion of the beam, are the negatives of F and 

n 

M found above. Fig. A.2.b. 
n 

PRINT -OUTS 

The print-outs are several pages as follows 
1st page. 

a. Pipe dimensions 

b. Computed cross sectional area, moment of inertia 
and polar moment of inertia of the cross section. 

c. Material properties, 

d. Number of points analysed in each cross section, 
number of cross sections analysed and their distances 
from the tip. 

2nd page. 

a. Radial distance of the points analysed . 

b. Loads applied to the tip. 

c. Distance from the tip to the first cross section 
analysed. 
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d. 



Computed loads acting at that cross section 



Following pages. For each analysed point the print-out contains: 

a. Angular position of the point in the cross section 
(degrees) 

b. Forces and moments acting at that cross section in 
radial, tangential and axial components. 

c. Stresses produced by individual loads. 

SHTZT is the shear stress in the ZT direction 
produced by the axial component of the moment 
SBZZ axial (ZZ direction) stress produced by 
the tangential component of the moment. 

SAZZ axial stress produced by the axial compon- 
ent of the force. 

SHFZR and SHFZT are stresses produced by the 
radial and tangential components of the force. 
SPRR and SPTT are radial and tangential stresses 
produced by the internal pressure. 

d. If the analysed points are in the external surface 
of the pipe, there will be additional print-outs of 
the total axial, tangential and shear stresses at 
the point. 

Analysis of other cross sections for the initial and/or other radial 
distance and tip loads and pressure produce similar outputs as described 
frcro page 2 on. The complete program and sample outputs of the analysis 
made for points at the external surface and % in. under that surface of 
cross sections 1 and 2 of the pipe shown in Fig. A.I., follow this page. 
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EXTERNAL diameter = 
INTERNAL DIAMETER = 
LENGTH 



CRUSSFCT I ONAL AREA = 
MOMENT JF INERTIA = 
POLAR MOMENT OF l . = 



MODULUS OF ELASTICITY 
SHEAR MODULUS OF ELAST 
COEFFICIENT OF THERMAL 
POISSON RATIO 0.300 
SPECIFIC WFIGHT 0. 



12.00 INCHES 

10.00 INCHES 
90. OC INCHES 



3 A . 5 5 7 5 I N 2 
527.0020 IN4 
1054.0039 IN4 



0.30E 03 PSI 
ICI TV O. 12F 08 PS I 

EXPANSION NOT GIVEN 

2 83 C LB/IN3 



NUMBER OF POINTS TO ANALYSE 12 

NUMBER OF SECTIONS TO ANALYZE 2 
DISTANCES F R CM THE TIP TO THF CROSS SECTIONS, 

30.000 60.000 0.0 0.0 C.O 



INCHES 
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RADIAL DISTANCE OP POINTS 



6.00000 INCHES 



LOADS APPLIED TO THE 

INTERNAL PRESSURE = 
APPLIED FORCE 
APPLIED MOMENT 



DISTANCE FROM THE TI 



LOADS ACTING ON THE 
INTERNAL PRESSURE = 
APPLIED FORCE 
APPL I ED ..MOMENT 



TIP 

100.00 PSI 

100.00 100.00 100.00 LB 

100.00 100. 00 100.00 LB-IN 



P TO THE CROSS SECTION 30. CO INCHES 



CROSS SECTION 

100.00 PSI 

393.39 100.00 100.00 LB 

-2900.00 7500.89 100.00 LB-IN 
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. " OF GF r ES 








T ►' <\SF. ) -• " 


> 


F - :CF = : .3*339= ' 3 0.100 TOP 


0 3 


P.lOOOPf C 3 LH 


1 s A f . 
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vpvu.xir =-C.29CC 


OF 09 0 . 7 5 0 o o p 
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ShI/7 
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— f> 


.357 S' 0 F 02 








5 A 7 7 
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SHI i' 




. P 
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s r :• •> 
5 L 7 T 


, > 


• Sj.'iMC Cl 

A 
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.1-99551 C 3 












7*7 P 1 r- F t 7 I UN 


Z'Z DIRECT ION 




Z'T DIRECTION 






AL l 1 ( 1 , 1,1) 


AIL 1(2, 1,1) 




ALLI ( 3, 1,1) 


S 7 F E 9 S 7 




~>.9 r w55E ' 3 


-0.R25O5F 02 




0 . 59 55 3E 01 


1 M~T \ = 


30 


. ' 3 DEGREES 








7 •-> !J 3 (- ■ > v 1 


j 


Ml- Ll = '.'.39069? 0 3-0. 1 1 0 A 9 F 
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7“ ,">F U’> 
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OF 09 P . 7997) J E 
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. 5 s 9 2 5 F rr 








S > 1 ■’ 
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7' 7 H I D F(, 7 I ON 


Z'Z DIRECT I f*N 




Z'T DIRECTION 






M L 1 ( 1 , 2,1) 


ALL 1(2, 2,1) 




A L L l ( 3 , 2,1) 


ST 1 1 S S c 3 




-7 . A 5 A 5 5 5 0 7 


-0.87572? 92 




-0.53605E 01 


r.i^TA = 


30 


. "0 OFGFFFS 








tr~nsf ,7-'.; 


0 


MlKCE = P.28330F 03-0.29089F 


03 


C . 10 00 OF 03 LB 


7 R A i 4 S F 7 \ ■ ' 1* 


) 


P i -/ iF NT = C . 50 4 6 9 F 09 0.62619? 


09 


0 . 190OOE 03 LB-IN 


SHT Z 7 


r 


. 5 6 9 2 <j F '.) C 








3 >71 




.71?9 7 F 02 








s tzz 


0 


.2393 7 F 01 








S H F 2 s 




A 

• 








SHT 7 T 




. 1 5 5 5 7 F C 2 








3 PR Is 


0 


.0 








S P 7 T 




. 9 5 9 5 5 F 03 












T ' T DIRECTION 


Z'Z 7IRFCTICN 




Z'T DIRECTION 






A L L 1 ( l , 7,1) 


ALLI ( ? , 3,1) 




ALLI ( 3 , 3,1) 


s r p r s s i_ s 




P.95955F 03 


- 0 . 53399? 02 




-0.150B7F 02 


7 Hi 7 A = 


9*: 


. DEGREES 








7 is - f S F u 7 v F 


r ) 


l PRC E = 0 . 10~( 


OF C 3-0. 3 9 33 9 F 


03 


P. 10009? 03 LB 


7 •> A 'i S F : !K M | 




H-HFNT = (.750C9E 09 0.?9000r 


09 
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SHT 7 T 
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. 3 3 0 1 7 F 0? 
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SH r Z T 
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3 PR 
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S P T T 


f' 


. 4-,4?SE 03 












T'T DIRECTION 


7. '7 DIRECTION 




Z'T DIRECTION 






AlLI ( 1 , 9,1) 


ALL 1( ? , 9,1) 




ALLI (3, 9, 1 ) 


T F h 3 3 c 3 




0.95955F 03 


-0.30123? 02 




-C.20M9F 02 
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THU A = 1 20 .00 DEGREES 

TRANSFORMED FORCE =-( . 11CC9E 03-0.39069F 03 C.10000E 03 L6 
T R A N S F 0 R y c 0 MOMENT = p . 794t>0E 04-0.12390E n 4 0.10000E 03 LR-IN 



SHTZT 
S nZZ 
S AZZ 
SHF ZR 
SHF Z T 

s r f •' 

S P T T 



0 . 5 < j 9 ? 6 E 00 
C.14106E 02 
C.23937E 01 
0.0 

-0.2l p 43E 0? 

0.0 

0.45 45 51- 03 



S T p F S $ E > 



T 1 T DIFFCTION 
ALL 1(1, 5,1) 

0.45455F 03 



1*7. DIRECTION 
ALL 1(2, 5,1) 

0. ) 6999E 02 



7 . * T DIRECTION 
ALL 1 ( 3 » 5,1) 
-0.20473E 02 



T H r 7 A - 150.00 DEGREES 

T '< A N SF OR ME o FORCE =-C.?«C69E 03-P.2833OE 03 0.10000E 03 LB 
TRAP S F • l IMP i' MOMENT = l.6261 p F C4-0.50463F 04 0.10000E 03 LB-IN 



SHTZ T 


P.5n926F 00 






S c 7 7. 


0.5744PF np 






S A 7 7. 


•'.2 39 3 7E 01 






SHF 79 


C . 0 






SHFZ T 


. 1 52 3 9E 0 2 
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n n 
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S p lT 


C.45455E 03 








T * T DIRECTION 


Z • 7 DIRECTION 


Z * T DIRECTION 




A L L 1 ( 1 , 4,1) 


AL 1 1(2, 6,1) 


ALL1 ( 3, 6,1 ) 


STRESSES 


0 . 45455F 03 


0.60343F 02 


-0.14689E 0? 



T HF T A = 


ISC.rj DFGRFES 




TRANSFORMED FORCE =-C.39339F 03~0. 10C00E 03 


C.10000F 03 L8 


TRANSFORM 


Ml! IF NT = F . 2 9 r 0 0 F 04-0.750OQF 04 


0 . 1 0000 E 03 LB-IN 


SHTZT 
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S rs Z 7 
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S A LI 


0.2 39j 7E 0 1 




SHF ZR 


p.n 




SHF ZT 


-0.5 386 IE r l 




S PRR 


c . c 




S PTT 


0.9 545 SE 03 






T * T DIRECTION 7*7. DIRECTION 


Z'T DIRECTION 




AL L 1 ( 1 , 7,1) A L l 1(2, 7,1) 


ALL1(3, 7,1) 


STP FSRES 


p .45455F 03 0.P8Z92F 0? 


-0.4B168E 01 



T H r T £ = 


210 


.0.' DFGPEFS 








Th AMSFrRV) 


0 0 


FORCF =-C.39069F 


03 


0.11039F 03 


0.10000F 03 LB 


T° ANSEtKMFO 


FOMENT =-0 . 12 390E 
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SHTZ T 
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S PTT 
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T • T DIRECT LUN Z 


• 7 


DIRECTION 


Z'T DIRECTION 






ALL 1 ( 1 , 0,1) 


All 1(2, 8,1) 


ALL1I3, 8,1) 

0.64989E 01 


STF [ SSES 




0.45455E 03 


0 . 


93 360 E 02 
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THFT \ = PA''. no DFGFFFS 

TRANSFORMED FORCE =-C . 28330E 01 0.29069F 03 O. 10000F 03 LR 
TRANSFORMED MOMENT =-C . 50 A ^9 E 0A-0.62619E CA 0.10000E 03 LB-IN 



SHT Z T 
S 07 7. 
S A Z Z 
SHF ZR 
SHF Z T 
S PRO 
S F'TT 



C.R6926F CO 
G.71293E 0? 
C.23037F r 1 
0.0 

p. 15657F r '2 

0.0 

p . A 5 A 5 5 E 03 



T • T DIRECTION 

ALL 1 II, 9,1) 

STRESSES O.A5A55E 03 



Z'Z DIRECT ION 
ALL 1(2, 9,1) 
0.7A187F 02 



7 . ' T DIRECTION 
ALL 1 ( 3 , 9,1 ) 
0.16226E 02 



THETA = 27c.O) DEGREES 

TRANSFORMED FORGE =-C. 1000CF 03 0.39339E 03 0.100C0E 03 LQ 
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0 . A5A55E 03 -0. 


5A555E 02 


0. 15828E 
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RADIAL DISTANCE OF POINTS 



6.00000 INCHFS 



LOADS APPLIED TO THE TIP 

INTERNAL PRESSURE = 
APPLIED FORCE 
APPLIED MOMENT 



100.00 PSI 

100.00 100.00 

100.00 100.00 



100.00 LB 
100.00 LB-IN 



DISTANCE FROM THE TIP TO THE CROSS SECTION 60.00 INCHES 



LOADS ACTING ON THE CROSS SECTION 
INTFRNAL PRESSURE = 100.00 PSI 

APPLIED FORCE * 686.79 100.00 

APPLIED MOMENT = -5900.00 23703.59 



100.00 L8 
100.00 LB-IN 
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THFTA = 


0.0 DEGREES 








TRANSFORMED 


FORCE = C.68679E 


03 


0.10000E 03 


C.10000F 03 LB 


TRANSFORMED 


MOMENT = -C • 59000E 


04 


0.23704E 05 


0.10000E 03 LB-IN 
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ALL 1 ( 3 , 4,2) 
-0 . 3642 IE 02 


STF FSSES 


0.45455F 03 -0.64279E 02 
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RADIAL DISTANCE OF POINTS 



5.50000 INCHES 



LOADS APPLIED TO THE TIP 

INTERNAL PRESSURE = 100.00 PS I 

APPLIED FORCE = 100.00 100.00 100.00 LB 

APPLIED MOMENT = 100.00 100.00 100.00 LB-IN 



DISTANCE FROM THE TIP TO THE CROSS SECTION 30.00 INCHES 



LOADS ACTING ON THE CROSS SECTION 
INTERNAL PRESSURE = 100.00 PS I 

APPLIED FORCE = 393.39 100.00 100.00 LB 

APPLIED .MOMENT- = -290Q..00 _ 7500.89 100.00 LB-IN 



TM! T A = 1 OKJttRES 

Tk/.KSFO^Mt') hOkCC =-r. llOC^E C3- 30069r 03 '-.lOOOCF 0 3 13 
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RADIAL DISTANCE OF POINTS 



5.50000 INCHES 



LOADS APPLIED TO THE TIP 

INTERNAL PRESSURE = 
APPLIED FORCE 
APPLIED MOMENT 



100.00 PSI 
100.00 
100.00 



100.00 

100.00 



100.00 LB 
100.00 LB-IN 



DISTANCE FROM THE TIP TO THE CROSS SECTION 



60.00 INCHES 



LOADS ACTING ON THE CROSS SECTION 
INTERNAL PRESSURE = 100.00 PSI 

APPLIED FORCE - 686.79 100.00 

APPLIED MOMENT ..... = ..-5900.00 .. 23703.59 



100.00 LB 
1 00.00 LB-IN 
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THf T \ 
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0. 


1 0 0 0 0 F 


03 


LB 


7 kANSF Vi ME 


) 


'•ii i v i r m t 


=-C .59roOE 


OA 


0. 237TAF 


05 


c . 


10000E 


03 


LB-IN 


SHT7 1 


r 


. 621H2F 


or 
















S h/7 


-0 


. 2 A 7 3 H F 


r> 3 
















S A7Z 


r. 


.20937F 


01 
















SHF 7 K 




.A601 iT f 


r, r 
















buezT 




. •> f 40?E 


0 1 
















S PPK 


- 


. a 320 1 F 


~2 
















s nr 


w 


. A 9 775E 


03 
















THF TA = 


3 ' 


. ''O DFGFFES 
















T PANSFfU -U 


) 


rORCE 


- r .64477c 


03- 


0. 26679F 


03 


c. 


1 O000F 


0 3 


t.B 


T9 AN S FH <*irO 


MOMENT 


= V.67A2PF 


04 


0. 2B473F 


A 5 


c. 


1000OF 


03 


LB-IN 


SHT 7T 


< ' 


. 62132E 


Or 
















S o 7 7. 


- 


. 2 4 6 0 ? E 


f 3 
















S A 7 7 


i ■ 


. 78°37F 


01 
















SHi- 7s. 


r 


. A2PSAR 


Of. 
















SF-f 7 r 


t 


. 1 A 7A T P 


0 ? 
















9 n k 9 


~ u 


.A HP'? IF 


H? 
















S PTT 




. AM77SF 


np 
















THfTA = 


6 l. 


.77 D 


t GPEES 
















TRANSFORMED 


F DkC c 


= G.4300GF 


03- 


O. SAA 7 7F 


03 


0. 


1 000 OF 


03 


1 B 


1 ft ...VSFi) -HFj 


Hfl-IEMT 


= f . 17578E 


06 


0. 1 6 0 6 1 E 


os 


0. 


10000F 


13 


LB- IN 


SMT 7 T 


p 


.^P18PF 


OC 
















S : 7 7 


- c 


. L 77CPF 


r 3 
















F A 7 Z 


c 


. 2 8°3 7F 


o 1 
















ShF 7P 


r* 


. 231P6F 


r< c 

V. % 
















SHF ZT 


-( 


. 31271F 


02 
















S PRO 


— i.* 


.4320 IF 


0? 
















S PIT 


( 


. 4 a 7 7 6 E 


o? 
















THfTA - 




. r O OfcGfttES 
















T R u \ S h j •s ^ i 


t ; 


t URGE 


= c.nccoe 


03- 


0.68679F 


03 


c . 


1 0 0 0 0 E 


0 3 


LB 


Tri ANSFORVF 


J 


Ml IMF 1ST 


= C .2370AE 


06 


0.B9Q00E 


04 


0. 


1 0000 E 


33 


LB- IN 


SHT7T 


l 


. c 2 1 H ? f 


,r 
















S o LL 


- 


16 7 3 t 


r, 7 
















S A 7 7 




. pnnp 7r 


01 
















SHF 7 V 




. a 5 5 A v< F 


- r 1 
















SHF 7 T 


- " 


. 3 9 A 2 3 F 


op 
















S *’F R 


- >' 


.A } ? 0 1 F 


' 2 
















c 'IT 




. - J 776F 


np 

















63 



THf^TA = 120.00 DEGREES 

TRANSFORMED FORCE =-0.25b79F 03-0.64477E 03 0.10C00E 03 IS 
TRANSFORMED MOMENT = 0.23478E 05-0.67422E 04 0.10000E 03 1.8-IN 



SHTZT 0.52182E 00 
S till C.70364E 02 
S A Z Z C.2SQ37E 01 
SHF ? R -0. 16832E OC 
SHE ? T - 0 . 3 701 IE 02 
S PER -C.A3201E 02 
S PTT C.49775F 03 



THETA = 15C.00 DEGREES 

TRANSFORMED FORCE =-C.54477E C3-0.43030F 03 0.10000E 03 L3 
TRANSFORMED MOMENT = 0. 16961E 05-0. 1 7573F 05 0. 10000F 03 L8-1N 



SHT7T 0.52182F OC 
S till C.13345F 03 
S A H C . 2 8 9 3 7 E 01 
SHE ZR -0.35709E OC 
SH T 7 T -0.24682F ^2 
S PPR -C.43201E 02 
S PTT C.49775E 0? 



THETA = 13C.00 DEGREES 

TRANSFORMED FORCE =— C. 68679F 03-0.10000F 03 0. 1000CE 03 10 
TRaNSFOTMFD MOMENT = C.590C0F 04-D.237D4F 05 C.1000CF 03 LB-IN 



SHTZT 0.521H2E 00 
S b 77 0.24738E 03 
S A 7Z C.28937E 01 
SHEZR -C.45C18E OC 
SHF 7 T -C . 57402 E 01 
S PR* -C. 4320 IF C? 
S PTT C.49775F 03 



THETA = 210.00 D E G REFS 

TRANSFIXED FORCE =-C . 644 77E C3 0.2F679F 03 C.10000E 03 LR 
TRANSFORMED MOMENT =-f.o7422E 04-0.23473E 05 C.10000F ">3 1.8- I N 

SHTZT F.52182F OC' 

S till C.24 5C2E 0.3 

S aZZ ' .28937E 01 

SHEZR -C.42264F OC 

SHTZT C.1474CE 02 

S PRR -C.432C1E 02 

S PTT C.A9775F 03 
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THf-TA = Z4C.n'J DEGREES 

TRANSFORMED FORCE =-0.4300GF 03 0.54477F 03 0.10C00F 03 LB 
TRANSFORMED Mi/M FN T =-f.L7578E 05-0.1*9616 ^5 0. 10000E 03 LB-IN 



SH17T C . 5 "> 1 8 2 E 00 
S dll 0. 177026 03 
S A 7 7 C.28937E 01 
SHF ZR -C.ZS186E OC 
SHFZT 0.3 127 IF 02 
S PRR -C.432C1F 02 
S PTT C. 497756 03 



THETA = 270. dd DFGREES 

TRANSFORMED F ° R C F =-C.100C0F 03 0.686796 03 0.10000E 03 LB 
TRANSFORMED MOMENT =-C. 237C4E 05-0.590^06 04 0.10000E 03 LB-IN 



SHTZT C.52192E OC 
S BZ 7 C.61575F 02 
S A Z Z C.2B937F 01 
SHF ZR -0.65 54 U F — 0 1 
SHF 7 T 0. 39423F 02 
S PRR -0.4 320 IF 02 
S PTT C.49775E ^3 



THFTA = 300.00 DEGREES 

TRmNSFQRMFO FORCE = 0.25679E 03 03 0.10000E 03 LB 

TRANSFORMED MOMENT =-C.23478E 05 0.67422E 04 O.IOCOOF 03 LB-IN 



SHTZT 0.52L82E OC 
S dll -0.70364E 02 
S A Z Z 0.2B937E 01 
SHF ZR C. 16832F OC 
SHFZT C.37011E 02 
S PRR -0.432C1F 02 
S PTT C.40775E 03 



THETA = 3 3t .00 DEGREES 

TRANSFORMED FORCE = 0.54477E 03 0.43000E 03 C.1000CE 03 LB 
TRANSFORMED MOMENT =-C.16961F 05 O.17570F 05 O.IOOOOE 03 LB-IN 

SHTZT 0.52182E CC 

S BZ Z -C.1B345E 03 

S AZZ C.2B937E 01 

SHFZR D.35709F 00 

SHFZT . 246826 02 

S PRR -C .4 3 2 0 1 E 02 

S PTT (. . 4 9 7 7 5 E 03 
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APPENDIX B 



DIGITAL COMPUTER PROGRAM FOR SECTION 2 AND NUMERICAL RESULTS 

The programs presented herein provide a method for applying the 
theory developed in Section 2 of this thesis. Actual situations can be 
handled by subroutines ZULIA, GUAYRA, CUMANA and DSIMQ (DSIMQ is present- 
ed in Appendix D) . Simulated situations require in addition subroutines 
CANTIL presented in Appendix A, and TIUNA and BLANCA presented in this 
Appendix. The last of the subroutines mentioned is an adaptation of a 
subroutine taken from reference 7 of this thesis. 

ZULIA is used to specify the required information, such as material 
properties, number of Fourier coefficients to be computed and descrip- 
tion of the arrangement of strain gage elements and their transverse 
sensitivity factor. Enough self explanatory comment statements have been 
included which specify the introduction of the information. Print-out 
statements present the information received for processing. Once the in- 
formation is properly arranged, subroutine GUAYRA is used to make the re- 
quired computations. 

The sequence in GUAYRA follows closely the theoretical development 
of the method. It has sufficient comment statements separating blocks 
of computer processing instructions as to follow the theoretical pro- 
cedure. Print-out statements present the computed Fourier coefficients. 

CUMANA is used by GUAYRA for auxiliary computations. 

Use of the subroutines is made by a very brief main program TERESITA. 
This program merely dimensions the various subroutines and turns control 
over to subroutine ZULIA. The card arrangement of the source deck is as 
shown in Fig. B. 1. 

Some sample results of the application of the program as specified 
at the end of Section 2, follow the card listings of the programs. 
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//TAB EXEC FORTCLG 
//TERESITA JOB '0324 JESUS A. TABORDA ' ,MSGLEVEL=1 



Fig. B.l 
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CALL ZULIA 
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THETA = 


0 


DFGREFS 


POINT 


STRESSES 




TANG. 




ACTUAL 




0.45455E 


C3 


COMPUTED 




0.A5A55E 


C3 


THETA = 


10 


DEGREES 


POINT 


STRESSES 




TANG. 




ACTUAL 




0.45455E 


C3 


COMPUTED 




0 .A5A55E 


C3 


THETA = 


20 


DEGREES 


POINT 


STRESSFS 




TANG. 




ACTUAL 




0. 4 545 5E 


C3 


COMPUTED 




0. A5A55E 


C3 


THETA = 


30 


DFGRFFS 


POINT 


STRESSES 




TANG. 




ACTUAL 




0. A5A55E 


C 3 


COMPUTED 




O.A5A5AE 


C3 


THETA = 


AO 


DEGREES 


POINT 


STRESSES 




TANG. 




ACTUAL 




0.45455E 


C3 


COMPUTED 




O.A5A5AE 


C3 


THETA = 


50 


CEGREES 


POINT 


STRESSES 




TANG. 




ACTUAL 




0.45455E 


C 3 


COMPUTED 




0 • 4 54 5 4 E 


C 3 


THETA = 


60 


DEGREES 


POINT 


STRESSES 




TANG. 




ACTUAL 




0. A5A55E 


C3 


COMPUTED 




0.A5A55E 


C3 


THETA = 


70 


DEGRFFS 


POINT 


STRESSES 




TANG. 




ACTUAL 




0. 45455E 


C3 


COMPUTED 




0. A5A55E 


C3 


THETA = 


80 


DEGREES 


POINT 


STRESSFS 




TANG. 




ACTUAL 




0 . A5A55E 


C3 


CGMPUTED 




0.A5A55E 


C3 


THETA = 


90 


DFGREFS 


POINT 


STRESSES 




TANG. 




ACTUAL 




0. A5A55E 


C3 


COMPUTED 




0. A5A55E 


C3 


THETA = 


100 


CEGREES 


POINT 


STRESSFS 




TANG. 




ACTUAL 




0. 45455E 


C3 


COMPUTED 




0 ♦ A 5A 5 5 E 


C3 


THETA = 


no 


CEGREES 


POINT 


STRESSES 




TANG. 




ACTUAL 




0. A5A55E 


C3 


COMPUTED 




0.A5A55E 


C3 



NO. 1 



LONG. 




SHEAR 




0. 17552E 


01 


0. 595 53F 


01 


0. 1755AE 


01 


0.59553E 


01 


NO. 11 








LONG. 




SHEAR 




0. 19702E 


01 


0.A9382E 


01 


0. 1970AE 


01 


0.49382E 


01 


NO. 21 








LONG. 




SHEAR 




0. 22133E 


01 


0.3788AE 


01 


0. 22135E 


01 


0.378 84E 


01 


NO. 31 








LONG. 




SHEAR 




0.24770E 


01 


0.25A07E 


01 


0. 2A773F 


01 


0.25407E 


01 


NO. A 1 








LONG. 




SHEAR 




0.27534E 


01 


0. 123 3 IF 


01 


0.27538E 


01 


0. 12331E 


01 


NO. 51 








LONG. 




SHEAR 




0.303A1E 


01 


-0.9A6 09E- 


■01 


0. 303A5F 


01 


-0.94599E- 


-01 


NO. 61 








LONG. 




SHEAR 




0. 3310AF 


01 


-0.1A022E 


01 


0. 33110E 


01 


-0. 14022E 


01 


NO. 71 








LONG. 




SHEAR 




0.357A2E 


01 


-0.26498F 


01 


0. 357A8E 


01 


-0.26A98E 


01 


NO. 81 








LONG. 




SHEAR 




0. 38172E 


01 


-0.37997E 


01 


0.38178E 


01 


-0 • 379 97 E 


01 


NO. 91 








LONG. 




SHEAR 




0.A0322E 


01 


-0.A8168E 


01 


0. A0328E 


01 


-0.48168E 


01 


NO. 101 








LONG. 




SHEAR 




0.A2126E 


01 


-0.567 02E 


01 


0. A2132E 


01 


-0.56702E 


01 


NO. Ill 








LONG. 




SHEAR 




0. A3530F 


01 


-0.6 33 41 E 


01 


0. A3535E 


01 


-0.633 4 1 F 


01 
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THETA = 120 
STRESSES 
ACTUAL 
COMPUTED 



DEGREES 

TANG. 
0.45455E C 3 
0.45455E C3 



POINT NO. 121 

LONG. 

0.44490E 01 
0.44495E 01 



SHEAR 

•0.67882E 

•0.67882E 



01 

01 



THETA = 130 
STRESSES 
ACTUAL 
COMPUTED 



DEGREES 

TANG. 

0.45455E C3 
0.45455E C3 



POINT NO. 131 

LONG. 

0.44977E 01 
0.44982E 01 



SHEAR 
•0.701 88E 
•0.701 88E 



01 

01 



THETA = 140 
STRESSES 
ACTUAL 
COMPUTED 



DEGREES 

TANG. 

0.45455E C3 
0.45455E C3 



POINT NO. 141 

LONG. 

0.44977E 01 
0.44983E 01 



SHEAR 
-0.701 88E 
-0.701 88 E 



01 

01 



THETA = 150 
STRESSES 
ACTUAL 
COMPUTED 



CEGREES 

TANG. 

0.45455E C3 
0.45455E C3 



POINT NO. 151 

LONG. 

0.44490E 01 
0.44496E 01 



SHEAR 

-0.67882E 01 
-0.678 8 2 E 01 



THETA = 160 
STRESSES 
ACTUAL 
COMPUTED 



DEGREES 

TANG. 
0.45455E C 3 
0.45455E C3 



POINT NO. 161 

LONG. 

0.43530E 01 
0.43537E 01 



SHEAR 

-0. 633 41 E 01 
■0. 633 41 E 01 



THETA = 170 
STRESSES 
ACTUAL 
COMPUTED 



CEGREES 

TANG. 

0.45455E C3 
0.45455E C3 



POINT NO. 171 

LONG. 

0.42126E 01 
0.42133E 01 



SHEAR 

■0.56703E 01 
-0.56702E 01 



THETA = 180 
STRESSES 
ACTUAL 
COMPUTED 



DEGREES 

TANG. 

0.45455E 

0.45455E 



POINT NO. 181 

LONG. 

C3 0.40322E 01 

C3 0.40328E 01 



SHEAR 
-0.481 68 E 
-0.48168E 



01 

01 



THETA = 190 
STRESSES 
ACTUAL 
COMPUTED 



CEGREES 

TANG. 

0.45455E C3 
0.454556 C3 



POINT NO. 191 

LONG. 

0.39172E 01 
0.38178E 01 



SHEAR 

•0.37997E 01 
-0.37997E 01 



THETA = 200 
STRESSES 
ACTUAL 
COMPUTED 



CEGREES 

TANG. 
0.45455E C3 
0.45455E C3 



POINT NO. 201 

LONG. 

0.35742E 01 
0.35747E 01 



SHEAR 

-0.26498E 01 
•0.26498E 01 



THETA = 210 
STRESSES 
ACTUAL 
COMPUTED 



CEGREES 

TANG. 

C.45455E 03 
0.45455E C3 



POINT NO. 211 

LONG. 

0.33105E 01 
0.33110E 01 



SHEAR 

•0.14022E 

•0.14022E 



01 

01 



THETA = 220 
STRESSES 
ACTUAL 
COMPUTED 



CEGREES 

TANG. 
0. 45455E 
0.45455E 



POINT NO. 221 

LONG. 

C3 0.30341E 01 

C3 0.30347E 01 



SHEAR 

-0.946 23E-01 
-0.946 18E-01 



THETA = 230 


CEGREES 


POINT 


NO. 231 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0.27534E 


01 


0. 12331E 


01 


COMPUTED 


0.45455E 


C3 


0.27541E 


01 


0. 1 23 3 1 E 


01 
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THETA = 240 


DEGREES 


POINT 


NO. 241 








STRESSES 


TANG. 




LONG. 




SH EAR 




ACTUAL 


0.454 5 5E 


C3 


0. 24770E 


01 


0. 25407E 


01 


COMPUTED 


0 . 45455E 


C3 


0.24777E 


01 


0.25407E 


01 


THETA = 250 


DEGREES 


POINT 


NO. 251 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0. 45455E 


C3 


0. 22133E 


01 


0. 378 83E 


01 


COMPUTED 


0.45455F 


C 3 


0.22139E 


01 


0. 378 83E 


01 


THETA = 260 


CEGREES 


POINT 


NO. 261 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 19702F 


01 


0 • 493 8 2 E 


01 


COMPUTED 


0.45455E 


C3 


0. 19708E 


01 


0.49382E 


01 


THETA = 270 


DEGREES 


POINT 


NO. 271 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 17552E 


01 


0.59553E 


01 


COMPUTED 


C.45455E 


C 3 


0. 17558E 


01 


0.5Q553E 


01 


THETA = 280 


DEGREES 


POINT 


NO. 281 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0. 45455E 


C 3 


0. 15748E 


01 


0 . 6 80 8 8 F 


D1 


COMPUTED 


0.45455E 


C3 


0. 15754E 


01 


0.68087E 


01 


THETA = 290 


CEGREES 


POINT 


NO. 291 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C 3 


0. 14345E 


01 


0.74726E 


01 


COMPUTED 


0.45455E 


C 3 


0. 14351E 


01 


0.74726E 


01 


THETA = SCO 


DEGREES 


POINT 


NO. 301 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0 . 4 54 55 E 


C 3 


0. 13385E 


01 


3. 792 67 E 


01 


COMPUTED 


0.45455E 


C3 


0. 13392E 


01 


0.79267E 


01 


THETA = 310 


DEGREES 


POINT 


NO. 311 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 12898E 


01 


0.815 73E 


01 


COMPUTED 


0.45455E 


C3 


0. 12905E 


01 


0.815 73E 


01 


THETA = 320 


CEGREES 


POINT 


NO. 321 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.454 55E 


C 3 


0. 12897E 


01 


0.81573E 


01 


COMPUTED 


0.45455E 


C3 


0. 12904E 


01 


0.815 73E 


01 


THETA = 330 


DEGREES 


POINT 


NO. 331 








STRESSES 


TANG. 




L ONG. 




SH EAR 




ACTUAL 


0.45455E 


C3 


0. 13385E 


01 


0. 79267E 


01 


COMPUTED 


0.45455E 


C3 


0. 13390E 


01 


0.79267E 


01 


THETA = 340 


DEGREES 


POINT 


NO. 341 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0. 45455E 


C 3 


0. 14345E 


01 


0 . 747 26 E 


01 


COMPUTED 


0.4 54 5 5 E 


C3 


0. 14348E 


01 


0 . 747 26 E 


01 


THETA = 350 


CEGREES 


POINT 


NO. 351 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 15748F 


01 


0.6808RE 


01 


COMPUTED 


0.45455E 


C3 


0. 15 75 IE 


01 


0.68088E 


01 
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KS IS AN INDEX INDICATING TYPE OF SOLUTION: 0 FOR NORMAL, I FOR SINGULAR 
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THETA = 


0 


DEGREES 


STRESSES 




TANG . 


ACTUAL 




0 .45455E 


COMPUTED 




0.45455F 


T H I TA = 


T ( 


DEGREES 


STRESSES 




TANG . 


ACTUAL 




0 . P 5 4 5 5 f 


COMPUTED 




0 .854551 


THtrTA = 


20 


DEGREE S 


STRESSES 




TANG. 


ACTUAL 




0.45455E 


COMPUTED 




0. 45455E 


THETA = 


30 


DEGREES 


STRFSSES 




TANG. 


ACTUAL 




0 .454551 


C 0 M E U T r: i ) 




0.45455F 


THETA = 


4 0 


(H GREFS 


STh ESSE^> 




TANG. 


ACTUAL 




0.45455E 


COMPUTED 




0.45455F 


THETA = 


50 


DEGREES 


STRESSES 




TANG. 


ACTUAL 




0.45455E 


COMPUTED 




0.45455E 


THETA •= 


6C 


DEGREES 


STRESSES 




TANG. 


ACTUAL 




0.45455E 


COMPUTED 




0.45485F 


IHETA = 


70 


CEGRFES 


STRESSES 




TANG. 


ACTUAL 




G.4 t >45 t 'E 


COMPUTED 




0.45455F 


THETA = 


80 


DEbkEFS 


STE E SSES 




TANG. 


ACTUAL 




0.45455E 


COMPUTED 




0 . 45455F 


THETa = 


9C 


DEGREES 


STRESSES 




TANG. 


ACTUAL 




0.45455E 


COMPUTED 




0.45455E 


THETA = 


IOC 


DEGREES 


STRESSES 




TANG. 


ACTUAL 




G .454 5 5 E 


COMPUTED 




0 .45455F 


THETA = 


1 10 


DEGREES 


STRESSES 




TANG. 


ACTUAL 




0.45455E 


COMPUTED 




0.454S6E 



POINT NO. I 

LONG. 

03 0 . I 7 552 E 

03 0.17556E 



POINT NO. 11 

LONG. 

03 0.19702E 

03 0.19706E 



POINT NO. 21 

LONG. 

03 0.22133F 

03 0.22137F 



POINT NO. 31 

LONG. 

G3 0.24770E 

03 0.24774F 



POINT NO. 41 

LONG. 

03 0.27534F 

03 0.27538E 



POINT NO. 51 

LONG. 

03 0.3G341E 

03 0.30345E 



POINT NO. 61 

LONG. 

03 0.33104E 

03 0.33109E 



POINT NO. 71 

LONG. 

03 0.36742E 

03 0.35746E 



POINT NO. PI 

LONG. 

03 0.38 l 72E 

03 0.38177E 



POINT NO. 91 

LONG. 

03 0.40322F 

03 0.40327E 



POINT NO. 101 

LONG. 

03 0.42126E 

C3 0.42132E 



POINT NO. Ill 

LONG. 

03 0.4353CE 

03 0.43535E 



SHFAR 

01 0.59553E 01 

01 0.59553E 01 



SHFAR 

Cl n . 49382 F 01 

01 0.49382E 01 



SHEAR 

01 0.37884E 01 

01 0.37884F 01 



SHEAR 

01 0.26407E 01 

01 0.25407E 01 



SHEAR 

01 0.12331E 01 

01 C.12331E 01 



SHFAR 

Cl -0 .94609 E— 0 1 

01 — 0 . 94604 E-Q 1 



SHEAR 

01 -0.14022E 01 

01 -0.14022E 01 



SHEAR 

01 -0.2649 8 E 01 

01 -0 . 26498 E 01 



SHFAR 

Cl -G.37997F 01 

01 -0.37997E 01 



SHEAR 

01 -0.481 68 E 01 

Cl -0.48168E 01 



SHFAR 

01 -0 . 56702 E 01 

01 -0.56702E 01 



SHFAR 

01 -0. 63341 E 01 

Cl -0.63341F Cl 
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THETA = 120 
STRESSES 
ACTUAL 
COMPUTET 



THETA = 130 

STRESSES 
ACTUAL 
COMPUTED 



THETA = 140 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 150 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 160 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 170 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 180 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 190 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 200 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 210 
STRESSES 
ACTUAL 
COM PUT ED 



THETA = 220 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 230 
STRESSFS 
ACTUAL 
COMPUTED 



DEGREES 

TANG. 

0.45455E 03 
0.45455E 0.3 



POINT NO. 121 

LONG. 

0 • 44490 E 01 
0.44495E 01 



SHEAR 

■0.67882E 01 
•C.67882E 01 



DEGREES 

TANG. 
0.45455E 03 
0 . 4 5455E 03 



POINT NO. 131 

LONG. 

0.44977E 01 
0.44983E 01 



SHEAR 

■0.70188E 01 
■0.70188F 01 



DEGREES 

TANG. 
0.45455E 03 
0.45455E 03 



POINT NO. 141 

LONG. 

0.44977E 01 
0.44983E 01 



SHEAR 

•0.70188E 01 
•0.70188E 01 



DEGREES 

TANG. 
0.45455F 03 
0. 45455 E 03 



POINT NO. 151 

LONG. 

0 . 44490 E 01 
0.44496E Cl 



SHEAR 

'0 . 67882E 01 
■0.67882E 01 



DEGREES 

TANG. 
0.46455E 03 
0.45455E 03 



POINT NO. 161 

LONG. 

0.4353CE 01 
0.43536E Cl 



SHEAR 

■0.63341 E 01 
■0.63341E 01 



DEGREES 

TANG. 
0.45455F 03 
0.45455F 03 



POINT NO. 171 

LONG. 

0.42126E 01 
0.42133E Cl 



SHEAR 

■0 . 567 03 E 01 
■0.56702E 01 



DEGREES 

TANG. 
0.45455E 03 
0.45455E 03 



POINT NO. 181 

LONG. 

0.4032 2E 01 
. 0.40329E Cl 



SHEAR 

■0.48168E 01 
■0.48168E 01 



DEGREES 

TANG. 
0.45455E 03 
0.45455E 03 



POINT NO. 191 

LONG. 

0.38172E 01 
0.38179E 01 



SHEAR 

-0.37997E 01 
■0.37997E 01 



OEGREES 

TANG. 
0.45455E 03 
0. 45455 E 03 



POINT NO. 201 

LONG. 

0.35742E 01 
0.35749E 01 



SHEAR 

■0.26498E 01 
■0.26498F 01 



DEGREES 

TANG. 
0.45455E 03 
0.45455E 03 



POINT NO. 211 

LONG. 

0.33105E 01 
0.33111F 01 



SHEAR 

•0.14022E 01 
■0.14022E 01 



DEGREES 

TANG. 
0.45455E 03 
0.45455E 03 



POINT NO. 221 

LONG. 

0.3C341E 01 
0.30347E 01 



SHEAR 

-0.94623E-01 

•0.94625E-01 



DEGREES 

TANG. 

0. 45455 E 03 
0.45455E 03 



POINT NO. 231 

LONG. 

0.27534E 01 
0.27541F 01 



SHEAR 

0.12331E 01 
0.12331E 01 
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THETA = 240 


l.)EGKt E S 


POINT 


NO. 241 






STRESSES 


TANG. 




LONG. 




SHEAR 


Ac 1 UAL 


0.45455E 


03 


0.24770E 


01 


0.2 5407 E 


COMPUTE.) 


0.45485E 


03 


G. 24777 E 


01 


0.25407E 


THETA = 2 SO 


DEGREES 


POINT 


NO. 251 






STRESSES 


TANG. 




LONG. 




SHEAR 


actual 


0.45455E 


03 


0.22133E 


01 


0. 37883E 


COMPUTE.) 


0.454S5L 


03 


0.22139E 


Cl 


0. 37883E 


IHETA = 26) 


OFbREES 


POINT 


NO. 261 






STRL SSES 


TANG. 




LONG. 




SHEAR 


ACTUAL 


0.45455F 


03 


0.1 9702E 


01 


0.49382E 


COMPUTE') 


0 .45455F 


03 


0. 19708E 


01 


0.49382F 


1 HE T •» = 2 70 


DEGREES 


POINT 


NO. 271 






ST R ESSE S 


T A N U . 




LONG. 




SHEAR 


AC TUAL 


0 . 4 54 5 5E 


03 


0. 17552E 


01 


0. 59553E 


COMPUTE.) 


0 .45455* 


03 


0. 17558E 


01 


0.59553E 


THETa = 280 


DFGREES 


POINT 


NO. 281 






STF ESSES 


TANG . 




LONG. 




SHEAR 


ACTUAL 


0.45455E 


03 


0. 15748E 


Cl 


0.68088E 


COMPUTE ) 


0 • 45455E 


03 


0. 15754G 


Cl 


0.68P87E 


ThtTA = 290 


DEGREES 


POINT 


NO. 291 
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EXPANSIGN WITH II FOURIER COEFFICIENTS 
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THETA = 24 0 


DEGREES 


POINT 


NO. 241 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0.24770E 


01 


0. 2 5407 E 


01 


COMPUTED 


0.44647E 


C3 


0.28761E 


01 


0. 275 06 E 


01 


THETA = 250 


DEGREES 


POINT 


NO. 251 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0. 45455E 


C3 


0.22133E 


01 


0.37883E 


01 


COMPUTED 


0.44979E 


C 3 


0.25514E 


01 


0 . 3 97 7 5 E 


01 


THETA = 260 


DEGREES 


POINT 


NO. 261 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 19702E 


01 


0.49382E 


01 


COMPUTED 


0.45222E 


03 


0.22239E 


01 


0.50755E 


01 


THETA = 270 


DEGREES 


POINT 


NO. 271 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 17552E 


01 


0 . 5 95 5 3 E 


01 


COMPUTED 


0.452SSE 


03 


0. 19561E 


01 


0.6C426E 


01 


THETA = 280 


OEGREES 


POINT 


NO. 281 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 15748E 


01 


0.68088E 


01 


COMPUTED 


0.45262E 


03 


0. 17848E 


01 


0.68838E 


01 


THETA = 290 


DEGREES 


POINT 


NO. 291 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 14345E 


01 


0.74726E 


01 


COMPUTED 


0.45243E 


C 3 


0. 17031E 


01 


0.75850E 


01 


THETA = 300 


DEGREES 


PCINT 


NO. 301 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 13385E 


01 


0.79267E 


01 


COMPUTED 


0.45347E 


C3 


0. 16666E 


01 


0.810 HE 


01 


THETA = 310 


DEGREES 


POINT 


NO. 311 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0.12898E 


01 


0.81573E 


01 


COMPUTED 


0.45576E 


C3 


0. 16207E 


01 


0.83706E 


01 


THETA = 320 


OEGREES 


PCINT 


NO. 321 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C 3 


0. 12897E 


01 


0.81573E 


01 


COMPUTED 


0.45836E 


C3 


0. 15346E 


01 


0.83472E 


01 


THETA = 330 


CEGREES 


PCINT 


NO. 331 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 13385E 


01 


0.79267E 


01 


COMPUTED 


0.459S3E 


C3 


0. 14234E 


01 


0.8C2 72E 


01 


THETA = 340 


OEGREES 


PCINT 


NO. 341 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 14345E 


01 


0 . 7 47 2 6 E 


01 


COMPUTED 


0.45968E 


03 


0. 13451E 


01 


0.74556E 


01 


THETA = 350 


DEGREES 


POINT 


NO. 351 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0. 45455E 


C3 


0. 15748E 


01 


0.68088E 


01 


COMPUTED 


0.45770E 


03 


0. 13746E 


01 


0.67035E 


01 
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KS IS AN INDEX INDICATING TYPE OF SOLUTION: 0 FOR NORMAL » I FOR SINGULAR 
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THETA = 240 


DEGREES 


POINT 


NO. 24 1 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0.24770E 


01 


0.25407E 


01 


COMPUTED 


0.44800F 


03 


0.27843E 


01 


0.27467E 


01 


THETA = 250 


DEGREES 


POINT 


NO. 251 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0.22133E 


01 


0.37883E 


01 


COMPUTED 


0.44886E 


03 


0.25321E 


01 


0.39935E 


01 


THETA = 260 


DEGREES 


POINT 


NO. 261 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


C.19702E 


01 


0.49382E 


01 


COMPUTED 


0.45C09E 


03 


0.22847E 


01 


0.51336E 


01 


THETA = 270 


DEGREES 


POINT 


NO. 271 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0. 17552E 


01 


0.59553E 


01 


COMPUTED 


0.45157E 


03 


0.20507E 


Cl 


0.61327E 


01 


THETA = 280 


DEGREES 


POINT 


NO. 281 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455F 


03 


0.15748E 


01 


0 . 6 8088 E 


01 


COMPUTED 


0.45314E 


C3 


0. 18395E 


01 


0.69618E 


01 


THETA = 290 


DEGREES 


POINT 


NO. 291 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0. 14345E 


01 


0.74726E 


01 


COMPUTED 


0.45465E 


03 


0. 16605E 


01 


0 * 75972 E 


01 


THETA = 300 


DEGREES 


POINT 


NO. 301 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0.13385E 


01 


0.79267E 


01 


COMPUTED 


0 . 455 92 E 


03 


0. 15229E 


01 


0.80215E 


01 


THETA = 310 


DFGREFS 


POINT 


NO. 311 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0. 12898E 


01 


0.81573E 


01 


COMPUTED 


C.45684E 


03 


0. 14343E 


01 


0.82239E 


01 


THETA = 320 


DEGREES 


POINT 


NO. 321 








STRESSFS 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


C3 


0. 12897E 


01 


0.81573E 


01 


COMPUTED 


0.45728E 


03 


0. 14007E 


01 


0. 82001E 


01 


THETA = 330 


DEGREES 


POINT 


NO. 331 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0.13385E 


01 


0.79267E 


01 


COMPUTED 


0.45722E 


03 


0.14253E 


01 


0.79524E 


01 


THETA = 340 


DEGREES 


POINT 


NO. 341 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0 . 45455E 


03 


0. 14345E 


01 


'0.74726E 


01 


COMPUTED 


0 .45666E 


03 


0. 15C90E 


01 


0.74891E 


01 


THETA = 350 


DEGREES 


POINT 


NO. 351 








STRESSES 


TANG. 




LONG. 




SHEAR 




ACTUAL 


0.45455E 


03 


0. 15748E 


01 


0.68088E 


01 


COMPUTED 


0.45567E 


03 


0. 16491E 


01 


0.68246E 


01 



100 



EXPANSION WITH 3 FOURIER COEFFICIENTS 
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THFTA = 24C 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 250 
STRESSES 
ACTUAL 
COMPUTED 



THFTA = 260 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 270 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 2 BO 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 29C 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 300 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 310 
STRESSES 
ACTUAL 
CGMPUTED 



THETA = 320 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 330 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 340 
STRESSES 
ACTUAL 
COMPUTED 



THETA = 350 
STRESSES 
ACTUAL 
COMPUTED 



DEGREES 

TANG. 

0.45455E 

0.45166F 



DEGREES 

TANG. 
0.45455E 
0. 4519PF 



DEGREES 

TANG. 
0.45455E 
0 . 4 52 3 1 E 



DEGREES 

TANG. 
0.45455E 
0. 45261 f 



DEGREES 

TANG. 

0.45455E 

0.45288E 



DEGREES 

TANG. 
0.45455F 
0 . 45 3 1 1 E 



OEGREES 

TANG. 

0.45455F 

0.45330E 



DEGREES 

TANG. 
0 .45455F 
0.45345E 



DEGREES 

TANG. 

0.45455E 

0.45353F 



DEGREES 

TANG. 

0.45455E 

0.45356F 



DEGREES 

TANG. 

0.45455E 

0.45353E 



OEGREES 

TANG. 

0.45455E 

0.45345E 



NO. 241 
LONG. 

0.24770E 01 
0.26709E 01 



NO. 251 
LONG. 

0.22133E 01 
0.24177E 01 



NO. 261 
LONG. 

0.19702E 01 
0. 21831E 01 



NO. 271 
LONG. 

0.17552E 01 
0.1 9741 E 01 



NO. 281 
LONG. 

0. 15748E 01 
0.17972E 01 



NO. 291 
LONG. 

0.14345E 01 
0.16577E 01 



NO. 301 
LONG. 

0.13385E 01 
0.15597E 01 



NO. 311 
LONG. 

0.12898E 01 
0.15064E 01 



NO. 321 
LONG. 

0.12897E 01 
0.14994E 01 



NO. 331 
LONG. 

0 . 1 338 5 E 01 
0.15388E 01 



NO. 341 
LONG. 

0 . 14345 E 01 
0.16234E 01 



NO. 351 
LONG. 

0.15748E 01 
0.17507E 01 



SHEAR 

0.25407E 01 
0.26847E 01 



SHEAR 

0.37883E 01 
0.39286E 01 



SHEAR 

0.49382E 01 
0.50736E 01 



SHEAR 

0.59553E 01 
0.60850E 01 



SHEAR 

0.68088E 01 
0.69321E 01 



SHEAR 

0.74726E 01 
0 . 75 890 E 01 



SHEAR 

0.79267E 01 
0.80358E 01 



SHEAR 

0.81573E 01 
0.82590E 01 



SHEAR 

0.81573E 01 
0.82518E 01 



SHEAR 

0.79267E 01 
0.80144E 01 



SHEAR 

0.74726E 01 
0.75540E 01 



SHEAR 

0.68088E 01 
0.68846E 01 



POINT 



POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 

POINT 

03 

03 
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APPENDIX C 



DIGITAL COMPUTER PROGRAM FOR SECTION 3 AND NUMERICAL RESULTS 

Subroutine MACUTO presented in this Appendix provides a method for 
applying the theory developed in Section 3. Data for simulated situa- 
tions can be obtained from subroutine CANTIL presented in Appendix A. 

MACUTO has been programmed following closely the development of the 
method. Self explanatory comment statements are included, which detail 
the type of data to be used and the procedural sequence of the method. 
Print-outs present the data processed and the computed loading. 

Use of the subroutines is made by a brief main program called 
TIBISAY. This program merely dimensions the subroutines, accepts and 
prepares the input data, and turns control over to subroutine MACUTO. 

The card arrangement of the source deck is similar to that shown in 
Fig. B.l. 

Some sample results of the application made of the programs as 
specified in Section 3 are presented following the card listings of the 
programs . 
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APPENDIX D 



ubroutine D3II.IQ of the Computer Facility Center of tho Naval Postgraduate 
chool, to solve linear simultaneous equations. 
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APPENDIX E 



Effects of Changing the Number of Fourier Coefficients and the 
Amount and Quality of the Data in the Application of the Theory Developed 
in Section 2. 

The theory developed in Section 2 of this thesis estimates Fourier 
coefficients of the distributions of axial, tangential, and shear stress- 
es assuming that enough data has been obtained by placing a sufficient 
number of strain gage elements around the external surface of a cross 
section of pipe* The correctness of the estimated coefficients depends 
on their number and on the amount and quality of the data. Depending on 
the relative orientation of the gage elements, they can pick information 
useful for determining the coefficients of all the distributions, or, if 
they are arranged in certain disadvantageous ways, they may be able to 
provide information for none* or only some of the coefficients. 

The requirement of inferring a given number N of Fourier coeffi- 
cients, demands having, at a minimum, a corresponding number M - M(N) 
of valid data items. If more data are available, there is a redundance 
of data which permits more accurate determination of the N coefficients. 

Specific rules for the necessary number and disposition of gage ele- 
ments to obtain a specified number of correct coefficients are not avail- 
able, but they probably could be worked out. Some idea of the gage ele- 
ments requirements can be obtained by testing actual cases and comparing 
results with correct ones obtained from other methods. The digital com- 
puter programs included in this thesis can be used as they are to perform 
such tests. The procedure tc follow Ls to simply specify in subroutine 
ZULIA the number of gage elements, their angular positions and angular 
orientations around the cross section, and the number of Fourier co- 
efficients to be computed. The print-out of this program presents the 
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distributions obtained in this way , shown adjacent to the theoretically 
correct ones, so that comparisons can be readily done. 

A preliminary survey of gage elements arrangements was made using 
that procedure, applied to data with no added random errors. Results 
obtained are presented following these pages. Many of these results 
might have been anticipated from theoretical considerations. The follow- 
ing nomenclature is used when applicable: 

(a) No. of coefficients 

3 to 11: the results presented were obtained using from 

three to eleven Fourier coefficients per series. These 
results apply to that number of coefficients. 

7: results apply only to seven coefficients. 

(b) Results 

good: indicates that the coefficients obtained describe 

correctly the tangential, axial and shear stress distri- 
butions . 

bad: The computed coefficients do not describe either 

of the distributions. 

only shear good: good coefficients were obtained only for 

the shear stress distribution. 

Inasmuch as the data were without added errors, the recovery should be 
and actually was excellent or no good at all; accordingly, there is 
nothing subjective in the judgment M good M or ,r bad'\ 
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Gage elements located at 12 points 
equally spaced 30° around the external 
circumference of the cross section. 



gage elements 
arrangement 
at each point 



No. of Fourier 

coefficients per series results 



45° 




T axis 



3 to 11 



good 




3 to 11 



good 




3 to 11 



only shear good 




3 to 11 



bad 




3 to 11 



bad 



3 to 11 



good 



3 to 11 



only shear good 




3 to 11 



bad 
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II) 



Gage elements located at 8 points 
equally spaced 45° around the external 
circumference of the cross section. 



gage elements 

arrangement No. of Fourier 

at each point coefficients per series results 




3 to 7 



good 



9 and 11 bad 

same 




3 to 7 



good 



same 



9 and 11 



bad 




3 to 7 



only shear good 



same 



9 and 11 



bad 




3 to 11 



bad 



ki 



bad 
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Ill) 



One gage element every 10° around the external 
circumference of the cross section, totaling 
36 elements. The first element at 0° from the 
tangential direction, and every consecutive 
one, incrementing this angle by 10°. 

from 3 to 11 coefficients; bad results 



One gage element every 10° around the external 
circumference of the cross section, totaling 
36 elements. The first element at 0° from the 
tangential direction, and every consecutive 
one, incrementing this angle by 50°. 

from 3 to 9 coefficients; good results 
11 coefficients; bad results 
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